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Abstract 

We present certain classical continuum long wave-length limits of prototype integrable quan- 
tum spin chains, and define the corresponding construction of classical continuum Lax op- 
erators. We also provide two specific examples, i.e. the isotropic and anisotropic Heisenberg 
models. 
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1 Introduction 



Locally interacting discrete integrable spin chains have been the subject of much interest since 
they cropped up in string theory in the study of the AdS/CFT correspondence [1]. Their 
classical, long wavelength limit, provides a connection to continuous u-models describing 
particular dynamics of the string (references on this subject can be found in e.g. [2, 3]). 

Our motivation for this work is to develop a Hamiltonian approach different in its principle 
from the usual Lagrangian formulation of the long wavelength limit, in order to use in 
cases where the latter cannot be applied. In our approach we start from the Hamiltonian 
integrability formulation (quantum R-matrix and Lax matrix) guaranteeing a priori Liouville 
integrability of the classical continuous models. This is done through a Lax matrix-classical 
r-matrix formulation, provided that some consistency checks be made. On all known specific 
examples it will be checked that it yields the same results as the Lagrangian approach. It is 
indeed a key result that the Poisson structure is the same, in all cases when comparison is 
available, as the canonical structure derived from the long wavelength classical Lagrangian. 
This thereby validates the procedure and allows to use it in more general situations where 
the Lagrangian approach may not be used, in particular as a systematic way to build more 
general types of classical continuous integrable models by exploiting the richness of the 
algebraic approach. 

This contribution is based on [4] , where the interested reader can find all the details of the 
construction. 

2 The general setting 

In this section we outline the general procedure for obtaining a classical Lax formulation 
from the classical limits of the R and monodromy matrices. 

A quantum c-number non-dynamical i?-matrix obeys the quantum Yang-Baxter (YB) equa- 
tion [5] 



where the labels i — 1,2,3 may include dependence on a complex spectral parameter Aj. 
The auxiliary spaces are in this case loop-spaces Vi ® C(Ai), where Vi are (isomorphic) 
finite-dimensional vector spaces. 
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Assuming that it! admits an expansion ( "semiclassical" ) in positive power series of a param- 
eter (usually denoted h) as 

i?i2 = 1 ® 1 + hru + 0{h^) , (2.2) 

the first non-trivial term arising when we substitute this in (2.1) is of order two and yields 
the classical YB equation 

[ri2, ns] + [ri2, rss] + [ri3, ^23] = . (2.3) 

This is the canonically known "classical Yang-Baxter equation". It is not in general the 
sufficient associativity condition for a classical linear Poisson bracket, except when r is non- 
dynamical and skew-symmetric (see e.g. [6]). We shall hereafter limit ourselves to such 
situations. 

A quantum monodromy matrix T is generically built as a tensor product over "quantum 
spaces" and algebraic product over "auxiliary space" of representations of the YB algebra 
associated to R. Namely, one assumes a collection operators assembled in matrices Lu, 
acting on "quantum" Hilbert spaces labeled by i and encapsulated in a matrix "acting" on 
the auxiliary space Vi. For any quantum space q they obey the quadratic exchange algebra 
[7, 8, 9] 

R12 Liq L25 — L2q Liq Ri2 , (2.4) 

where operators acting on different quantum spaces commute. The form of the monodromy 
matrix T is then deduced from the co-module structure of the YB algebra 

Ta = Lai La2 ■ ■ ■ LaN (2-5) 

and thus naturally obeys the same quadratic exchange algebra (2.4). In particular one can 
pick L — R, the operators now acting on the second auxiliary space identified as "quantum 
space". This way, one builds closed inhomogeneous spin chains with general spins at each 
lattice site (labeled by (i)) belonging to locally chosen representations of some Lie algebra 
(labeled by i). 

Wc now establish that T has a classical limit by considering in addition the classical counter- 
part of L, labeled by L'^ which then satisfies the quadratic Poisson algebra, emerging directly 
as a semi-classical limit of (2.4), after setting ^[A, B] {A, B}. It reads 

{LUX,), LliX^)} = MX, - A2), LUX,) LliX^)] . (2.6) 
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The quantum monodromy matrix has also a classical limit given by (see also [10, 11]) 

Ta,{i} = ... I/^AT . (2.7) 

The exchange algebra for T'^ takes the form 

{T„^Ta = [ra6, . (2.8) 

This quadratic Poisson structure implies that the traces of powers of the monodromy matrix 
tr{T'^) generate Poisson-commuting quantities identified as classically integrable Hamiltoni- 
ans. Performing the trace over the finite vector space yields a generating function tr(T'^(A)) 
for classically integrable Hamiltonians obtained by series expansion in A. 



2.1 The long wavelength limit 

The usual presentation of the long wavelength limit, such as that found in [2, 3], is a La- 
grangian one where the Poisson structure is obtained from the standard derivation of canon- 
ical variables using a Lagrangian density. Instead, we will present here a purely Hamiltonian 
version of this limit by defining the long wavelength limit of a hierarchy of integrable quan- 
tum Hamiltonians based on some affine Lie algebra G. We shall define a priori the Poisson 
structure of the classical variables by imposing classical integrability of the long wavelength 
limit of the Hamiltonian through its associated classical Lax matrix. We consider a A'^-site 
closed spin chain Hamiltonian initially assumed to be governed by a nearest- neighbour 
interaction that takes the form 

N 

H = Y,Hii+i. (2.9) 
1=1 

The classical, long wavelength limit, is obtained by first defining local quantum states as 
linear combinations of the base quantum states. The bras and kets are denoted respectively 
by {n{l,9k)\ and \n{l,9k))-, where / denotes the site index and the dk^ denote the set of k 
angular variables. The condition of "closed" spin chain, essentially formulated as + / = /, 
imposes periodicity or quasi-periodicity conditions on the ^^^'s. Note that we have assumed 
that the base quantum states differ only by the fact that they are defined in distinct sites, 
hence the frequently used notation below instead of \n{l, 9k)), should not be confusing. 

If one considers nearest-neighbor local interactions then one defines the classical, but still 
defined on the lattice, Hamiltonian as 

N 

n = ^Uiit) , ni{x,t) = {ni\ ® Hu^i \ni) ® In^i) ■ (2.10) 

1=1 
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For integrable models, we may similarly define the continuum limit of the full set of commut- 
ing Hamiltonians. In these cases the generic Hamiltonians if^"^ of the integrable hierarchy 
are obtained directly from the analytic series expansion around some value Aq of the spectral 
parameter of the trace of the monodromy matrix (transfer matrix) as 

oo 

trT(A) = J](A - Ao)"//(") . (2.11) 

n=l 

By extension, we define in this case the classical Hamiltonians as the expectation value, over 
the N site lattice quantum state, of i?^") 

n^-\x,t) = (^f . . . {ml {ni+,\ . . . . . . \ni) \ni+i) .... (2.12) 

We next define a continuous limit and take simultaneously the thermodynamical limit in 
which iV — 7> oo. Accordingly, this is achieved by identifying the lattice spacing d as being 
of order 1/A^ and subsequently consider only slow- varying spin configurations (the long 
wavelength limit proper) for which 

k^l{x), k+i^l{x + 5). (2.13) 

In this limit, the finite "site differences" turn into derivatives. 

Given that (2.12) is applied to Hamiltonians of the integrable hierarchy obtained directly 
from the series expansion of the trace of the monodromy matrix, it is immediate that the 
expectation value procedure goes straightforwardly to the full monodromy matrix T (and 
thence to its trace over the auxiliary space which is altogether decoupled from the quantum 
expectation value procedure). Accordingly, we define first a lattice expectation value 

Ta = ...{nil® {Lai . . . LaN) ■ ■ ■ \ni) ® \ni+i) . . . , (2.14) 

which nicely factors out as 

N 

i=l 

Assuming now that L admits an expansion in powers of 5 as 

La^ = l + 6la^ + 0{6^) , (2.16) 

we consider the product (setting {ni\lai\ni) — la{xi)) 

N oo 

Ta = l[(^ + Slai + J2^-&- (2-17) 

i=l n=2 
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Expanding this expression in powers of 5, we get 

Ta = l + 6J2la^ + S'J2la,laJ+S'J2l'2 + ... ■ (2.18) 

i i<j i 

These, muhiple in general, infinite series of the products of local terms, are characterized 
by two indices: the overall power n of d, and the number m of the set of indices i (that 
is the number of distinct summation indices) over which the series is summed. Note that, 
in the T expansion one always has n ^ m. The continuum limit soon to be defined more 
precisely, will entail the hmit 5 — > with 0{N) — 0{l/5). We now formulate the following 
power-counting rule, that is terms of the form (for notational convenience lai — l^^^ below) 

m 

E J2^J-ri, (2.19) 

ii<i2<...im i=l 

with n > m are omitted in the continuum limit. For a sinlge summation, the latter is defined 

by 

Sy2la^^ / dxlaix) (2.20) 
i 

and similarly for multiple summations. Here A is the length of the continuous interval defined 
as the limit of A'^^. In other words, contributions to the continuum limit may only come 
from the terms with n = m for which the power 6"' can be exactly matched by the "scale" 
factor N"^ of the m-multiple sum over m indices i. In particular, only terms of order one in 
the S expansion of local classical matrices Lai = (^^il-^ail^i) will contribute to the continuum 
limit. Any other contribution acquires a scale factor — >• 0, when the continuum limit is 
taken. This argument is of course valid term by term in the double expansion. Being only 
a weak limit argument, it always has to be checked for consistency. 

Let's remark that if L is taken to be R, one naturally identifies S with the small parameter 
h, thus identifying in some sense the classical and the continuum limits. However, this is not 
required in general. It is clear to characterize separately both notions in our discussion as 

classical hmit : R— 1 + hr , 

continuum limit : L — 1 + 61 . (2.21) 

Recalling (2.13), the continuous limit of T, hereafter denoted T, is then immediately iden- 
tified from (2.15), as the path-ordered exponential from x to x = A 

r = Pexp (^J^ dx l{x)^ , (2.22) 
5 



where suitable (quasi) periodicity conditions on the continuous variables 9k{x) of the classical 
matrix l{x), acting on the auxiliary space V<S>C{X), are assumed. Of course the definition of 
a continuous limit requires that the L-matrices are not too inhomogeneous (e.g. L-matrices 
at neighbor sites should not be too different). This is in fact assured by the long wavelength 
limit assumption. 

2.2 The Lax matrix and r-matrix formulation 

The above identification of T also defines it as the monodromy matrix of the first order differ- 
ential operator d/dx + l{x). In addition, it has been built so as to straightforwardly generate 
the classical continuous limit of the Hamiltonians in (2.12) from the analytic expansion 

oo 

tr(r(A)) = ^(A - XoTn^""^ . (2.23) 

n=l 

We thus characterize l{x) as a local Lax matrix yielding the hierarchy of continuous Hamil- 
tonians 'H^'^\ In order for this statement to agree with the key assumption of preservation 
of integrability we are now lead to require a Poisson structure for I (inducing one for the 
continuous dynamical variables 9k{x)) compatible with the demand of classical integrabil- 
ity of the continuous Hamiltonians. Indeed, such a structure is deduced from (2.6), as the 
ultra-local Poisson bracket 

{h{x, Ai), h{y, A2)} = [ri2(Ai - A2), h{x, Ai) + h{y, A2)]5(a; - y) , (2.24) 

where r is the classical hmit (2.2) of the i?-matrix characterizing the exchange algebra of the 
L-operators. More specifically, recalling that Lai = 1 + Slai + 0{5'^), plugging it into (2.6) 
and assuming ultra-locality of Poisson brackets one gets 

{Li, hj} — [fab, Li + hj]^ ■ (2.25) 

One then identifies, in the continuum limit S ^ 0, the factor Sij/S with S{x — y). We 
then obtain a hierarchy of classically integrable, mutually Poisson commuting Hamiltonians 
from the explicit computation of the transfer matrix t{\) of the Lax operator d/dx -\- l{x) 
as if^") = ^^(A)|a=Ao- Such Hamiltonians are however generally highly non-local and not 
necessarily very relevant as physical models. We shall thus extend our discussion to local 
Hamiltonians. 
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2.3 Local spin chains 



Local spin chain Hamiltonians are more interesting, physically meaningful and easier to 
manipulate. In particular, they are the most relevant objects in connection with string theory 
and the AdS/CFT duahty [1]. Their construction generically requires the determination of 
a so-called "regular value" Aq of the spectral parameter such that Laj(Ao) oc Vai, where V 
is the permutation operator. In this sense the expansion of L can be expressed up to an 
appropriate normalization factor as 

L(A) = /(A)(l + 51 + 0(5')) . (2.26) 

Of course only when the auxiliary space a and quantum space i are isomorphic has this 
"regular value" any relevance. One then defines the local Hamiltonians as (denoting as 
usual t{\) = tr„r„(A)) 

H'-"' - ^Mm.^., ^ (2-27) 

Their long wavelength limit (e.g. (2.10)) is not obviously derivable from a straightforward 
"diagonal" expectation value of the T- matrix contrary to (2.12), since in general {F{A)) ^ 
F{{A)), for any functional of a set of operators A. However, we show below that this is 
indeed the case due to locality properties. Let us first focus for simplicity (but, as we shall 
see, without loss of generality) on the first local Hamiltonian 

^''' = '^^°^"^'(^^l^-Ao' (2-28) 

where, t^^(Ao) = 'Pi2'P23 • • ■'Pn-in- This operator acts exactly as a one-site shift on ten- 
sorized states, identifying of course site labels according to the assumed periodicity, i.e. 
N + 1 — 1. (Normalization issues are discussed in [4]). Computing the expectation value of 
we obtain 

, N 

= (ml (g) . . . (g) (n^|ri(Ao) — (/^(A)rr, + 5lai + 0{S^))'^ \n,) (g) . . . (g) M2.29) 

1=1 

One has 

(nil (g) (n2| . . . (njvK~"^(Ao) = (n2| {n^l (g) . . . (ni| (2.30) 
and of course iV -h 1 = 1. 
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Taking into account the power- counting rule described in section 2.2 we obtain that 



N r N 



(i/W) = n(n,+i|n,) J^(r(A)tr„n(l + 5( ) + 0{S'))) . (2.31) 

i=l 1=1 

We then easily cstabUsh that in the continuum hmit, using the power counting rule and the 
factorized form of both the state vector as ^ . . . ^ (n^l that the operator to be valued 
over it t'^iXo) = ^12^23 • • • ^iv-iiv, (^"H^)) = (^('^o))"^- We finally obtain that in the 
continuum limit 

{H^'') - r^(Ao)^t(A)|,^,; = (^(Ao))-^^(^(A))L^,„ = ^(in(t(A)))|,^,^ . (2.32) 

The computation may be easily generalized along the same lines for any higher Hamiltonian. 
Higher local Hamilltonians are indeed obtained from (2.27), admitting thus an expansion as 

= r^(Ao)^t(A)|,^ +polynomials , (2.33) 

depending only on lower order local Hamiltonians. When computing the expectation vahie 
of such higher Hamiltonians one gets the expectation value of t~^(Ao)^t(A)|Ao which in the 
continuum classical limit yields 

r'(Ao)^t(A)U„) = {tiXo))-'-^m)\x, , (2.34) 

using the same arguments as in the n = 1 case. In addition, one obtains expectation values 
of the polynomials of order k in the local Hamiltonians. In this case expectation values 
by tensor product of local vectors (ni| . . . {unI are exactly factorized over products of k 
local monomials hi^ . . . hk^ , except if some of the indices i coincide (or at least overlap for 
multiple indices). Locality of the lower Hamiltonians plays here a crucial role. It is clear that 
such families of terms with coinciding or overlapping indices correspond to a second "label" 
M — k — 1 and therefore their contribution will necessarily be suppressed in the continuum 
limit, with respect to the contribution of the generic terms (non-coinciding indices) with 
M = A; by the power-counting argument. Hence, it is consistent to conclude that in the 
continuum limit 

(Polynomial in {H^^)) = Polynomial in {{H^''^)) (2.35) 

and therefore 



(^^"^) = (:7T7rM^(A)). . ) = TTrln((^(A))) 



^dA" 



A=Ao dX^' 



A=Ao 



(2.36) 



This is the final, key result in systematically establishing the classical continuum limit of 
integrable spin chains. We may now apply this general procedure to all sorts of examples, 
starting with the simpler applications. 
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3 Examples 



3.1 The isotropic Heisenberg model 

The isotropic Heisenberg model (XXX chain) Hamiltonian describing first neighbor spin-spin 
interactions is given by 

^ = 2 E (-J^J+i + + -Mh-i) • (3-1) 

It is well known that when one considers the long wavelength limit one obtains a classical 
(j-model [2, 3]. We shall briefly review how this process works. The coherent spin state is 
parametrized by the parameters x, t via the fields 9, (p as 

\n{x,t)) = cose{x,t) e*'^^^'*) |+) + sin^(a;,t) e^^^^^^'*) |-) , (3.2) 

where the ranges of variables is ^ e (0, 7r/2) and (p e (0, tt). One can verify the completeness 
relation 

J dfx{n)\n){n\^l , (3.3) 

where the integration measure is given by 

4 

dfi{n) = — sin 6* cos 6' d9 dip . (3.4) 

TT 

Then as was described in [2, 3] and in subsection 2.1, one obtains a classical Hamiltonian 
via the expectation value procedure by employing (2.10). The appropriate XXX 2-site 
Hamiltonian is 

Hii+, oc {Vii+i - I) , (3.5) 

where V is the permutation operator acting as V{a (S>b) — b<S> a ior a, b vectors in V. Prom 
the definition of T-L we are led to compute quantities of the type 

(a| ^ {b\ V \a) O \b) = {a\b) ^ {b\a) = \{a\b)\^ . (3.6) 

They are expressed in terms of scalar products of the form 

(n|n) = cos(^ — 9) cos{ip — ip) + i cos{9 + 9) sva.[ip — (p) . (3.7) 
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In the long wavelength limit, \n) — \h) — \5n), 9{x) — 9{x + 6) and (p{x) — (p{x + 6). We 
finally conclude that 

Ho^ j dx {6'^ + sm\2e) ^'^) . (3.8) 

We shall now derive the Lax representation yielding (3.8) following section 2. The i?-matrix 
for the XXX model is [12] 

R(X) = x + inv . (3.9) 

This i?-matrix is a solution of the quantum YB equation [5] . It has a consistent normahzed 
classical limit defined as 

r(A) = \v , (3.10) 

which satisfies the classical YB equation. Alternatively, the classical r-matrix may be written 
as 

Set first 

^an(A) = i?„„(A-^) (3.12) 

and demand that L satisfies the fundamental algebraic relation 

Rabi^l ~ ^2) Lani^l) -^671(^2) = -^^(^2) -^an(Al) Rab{^l ~ -^2) , (3.13) 

where as usual in the spin chain framework we call n the quantum space and a the auxiliary 
space. Following the general derivation of section 2 and going directly to the continuous limit 
we disregard higher powers in S — h (in this case the two small parameters are naturally 
identified) . We next define a "local Lax matrix" as a mean value of L on the same coherent 
spin state, taken solely over the quantum space 

{n\Lan{\)\n) ^l + ihl{x,\) , (3.14) 

where 

'|(n|cT^|n) {n\a~\n) \ _ 1 / cos2^(a;) sin2^(x) e'^'^^^^ . 
{n\a~^\n) —^{n\a^\n) I 2 \ sin 29 (x) e+2^'^(^) -cos2^(a;) ' ' 
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where we have used the form of the coherent states to compute the matrix elements exphcitly. 
Then I satisfies the classical fundamental algebraic relation 

{li(x,Xi), Z2(z/,A2)} = [ri2(Ai-A2), li(Xi) + l2(X2Mx - y) . (3.16) 

Setting l{x,X) — 11/ A and taking into account the above algebraic relations we get 

{Hi, U2} = Vi2{ll2 - Ili)S{x - y) . (3.17) 

The parametrization in terms of the continuum parameters 9{x), (j){x) gives rise to the 
classical version of 5/2- Indeed, parametrizing the generators of the classical current algebra 

as 

= cos 29 , 2^ e^2^<^ . (3.18) 

we obtain from the fundamental relation that 

{S+,S-} ^ S'S{x-y) , {S',S^} ^±2S^S{x-y) . (3.19) 

The continuum parameters 6{x) and (j){x) can also be expressed in terms of canonical vari- 
ables p and q as 

cos29{x) = p{x) , (p{x) = q{x) and {q{x), p{y)} = i6{x — y) . (3.20) 

The Z-matrix in (3.15) coincides obviously with the potential term in the Lax matrix of the 
classical Heisenberg model. Precisely, one recalls that one must consider as classical Lax 

operator a la Zakharov-Shabat L = d/dx + l{x). The monodromy matrix for L is well 
known now to yield the classical Hamiltonians including the first non trivial one (see [10]) 

Recalling (3.18) and substituting in the expression above we obtain the Hamiltonian (3.8), 
hence the process above works consistently. 

Having exemplified the general construction of Section 2 to a simple system and checked the 
consistency of the approach we now turn to more complicated systems by first moving to 
trigonometric and elliptic sl{2) i?-matrices, corresponding to the XXZ and XYZ spin chains. 



(3.21) 
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3.2 The anisotropic Heisenberg model 



Consider the generic anisotropic XYZ model with Hamiltonian 



1 ^ 



For the following computations it is convenient to set 



(3.22) 



= 1 - S^a^ , C^{x, y, z} . 



(3.23) 



The Hamiltonian is written as 



N 



N 



y y 



N 



2 



(3.24) 



The additive constant may be omitted here. Taking into account equations (3.5)~(3.8), 
(3.24) and keeping terms of order S"^ we get 

H (X J dx {e"^ + sin2(2^) + sm^{2e) cos^(299) + ay siii^{2e) sin^(2(^) + cos^(2^))(3.25) 

This may be seen as an anisotropic "deformation" of the classical Heisenberg Hamiltonian. 
The last three terms are essentially potential- like terms. In the special case of the XXZ model 
the terms with coupling constant a^u o-y ai'e zero, whereas in the XXX case all potential terms 
vanish and one recovers the Hamiltonian (3.8). If we now recall the parametrization (3.18), 
then the expression above reduces to the Hamiltonian of the Landau-Lifshitz model or the 
anisotropic classical magnet [10] 



(3.26) 



We now derive the classical Z-matrix for the anisotropic cases. We focus in more detail on 
the XXZ i?-matrix 




R{\) 



sinh(A + f + f ) sinh(i/x)cT- 



sinh(i/x)(7"' 



sinh(A - f (7^ + f ) 



(3.27) 



The classical limit of the XXZ i?-matrix, after appropriate normalization, is given as (we 
divide with the constant factor sinh A) 



R{X) = 1 + i/x r(A) + 0{n^) , 



(3.28) 
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where 

r(A) = ^f^^ + ^^,^°^'' . \ \. (3.29) 

^' sinhA \^ (7+ (-^ + i) cosh Ay ^ ' 

The associated classical Lax operator is again obtained from L(A) = i?(A — ^) as (once again 
moving immediately to the continuous limit) 

{n\L{\)\n) + l(x, A) + 0{iJ^) , (3.30) 



with 



/(A) 



{n\^\n) cosh A {n\a \n) 
sinhA \^ (n|(7+|n) —{n\^\n) coshA^ 



1 (^S'coshX S- 



sinhA V S+ -^S' cosh A 



(3.31) 



where , are the classical generators of the current sl{2) algebra reahzed in terms of the 
angular variables in (3.18). The continuous variables x,y were omitted here for simplicity 
and will be from now on whenever there is no ambiguity. 

Let us also briefly characterize the classical algebra underlying the model. We set 

, cosh A ^ 1 , coshA^ 1 , , , 

smhA smh(A) smhA smh(A) 

Substituting this expressions to (3.16) and taking into account that 

[Ai2, A,] = -[Pi2, A2] , (3.33) 
we end up with the following set of Poisson structures 

{Di, D2} = Q, {Di, A2} = [D12, A2]6{x - y) , {Ai, A2} = -[A12, Di]6{x - y) (3.34) 
which give rise to the 5/2 Poisson algebra (3.19). 

The full XYZ classical r-matrix also yields, through this process, the classical Lax operator 
of the fully anisotropic classical Heisenberg model, satisfying also the fundamental linear 
algebraic relation (3.16) (see also [10]). A detailed presentation of this derivation is omitted 
here for the sake of brevity. 

Generahzations of the Heisenberg model associated to higher rank algebras are also presented 
in [4], where a more detailed analysis of the described process can be found. 
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